Dynamic Phases of Vortices in Superconductors with Periodic Pinning 
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We present results from extensive simulations of driven vortex lattices interacting with periodic 
arrays of pinning sites. Changing an applied driving force produces a rich variety of novel dynamical 
plastic flow phases which are very distinct from those observed in systems with random pinning 
arrays. Signatures of the transition between these different dynamical phases include sudden jumps 
in the current-voltage curves as well as marked changes in the vortex trajectories and vortex lattice 
order. Several dynamical phase diagrams are obtained as a function of commensurability, pinning 
strength, and spatial order of the pinning sites. 



PACS numbers: 74.60. Ge, 74.60. Jg 
The dynamics of driven vortex lattices interacting with 
quenched disorder have recently attracted considerable 
attention. Both theoretical and experimental ||] 

studies have suggested the exciting possibility that these 
systems exhibit novel dynamic phase transitions between 
different flow regimes as a function of driving force. In 
samples with random pinning, evidence for a transi- 
tion from a plastic flow phase to a strongly driven or- 
dered homogeneous phase has produced current debate 
over whether the strongly driven phase forms a moving 
crystal Q or moving ordered glass The dynamic 

phases may also be relevant to many other systems, such 
as Josephson junction arrays, charge-density waves, and 
electron crystals. Although much work has been done 
on dynamic phases in systems with random disorder, the 
case of periodic pinning has not been addressed. We show 
that driven vortex lattices interacting with periodic pin- 
ning exhibit a number of novel plastic flow phases which 
are not observed in random pinning arrays. Further, the 
onset of these different phases produces microscopic fea- 
tures in the vortex structure and flow patterns and gives 
rise to very pronounced features in macroscopically mea- 
surable current-voltage curves. 

General interest in periodic arrays of pinning sites has 
increased now that it is possible to construct samples 
with well defined periodic pinning structures in which 
the microscopic pinning parameters, such as size, depth, 
periodicity, and density, can be carefully controlled [j^ |) . 
Interesting commensurability effects are observed both in 
magnetization measurements Q and with direct imaging 
Q . Periodic pinning arrays are also of technological im- 
portance since they can produce higher critical currents 
than an equal number of randomly placed pins (9) . 

When the pinning radius is much smaller than the lat- 
tice spacing, the "disorder" in the system can be fine 
tuned by changing the commensurability. At B/B^ — 1, 
where B^ is the field at which the number of vortices N v 
equals the number of pinning sites Np, the vortex lattice 
locks into a periodic pinning array R~hI an d the pinning 
force is maximized. For B/B^ > 1, the vortex lattice 
contains two species of vortices: the pinned vortices that 



are commensurate with the pinning array, and the gen- 
erally more weakly pinned interstitial vortices that are 
caged by vortices at the pinning sites. When B/B^ < 1, 
the vortex lattice contains a well defined number of va- 
cancies which have their own dynamical behavior. 

In order to study the two-dimensional (2D) dynam- 
ics of rigid flux lines driven over periodic arrays of pin- 
ning sites, we have performed a large number of T = 0, 
current-driven molecular dynamics (MD) simulations. 
Unlike previous current-driven simulations (l],|^,^|,|]] , we 
examine the effects of periodic pinning arrays rather than 
random arrays, and cover a much larger range of the mi- 
croscopic pinning and system parameters, allowing us to 
construct a series of detailed dynamical phase diagrams. 

We numerically integrate the overdamped equations of 
motion (10): f j = f f v + f- p + t d = ijVi. Here f* is the 
total force acting on vortex i and we take rj = 1. The 
force from the other vortices is f™ = Ylj=i fo Ki(\ r i — 
Tj|/A) fy, where K\(r/X) is a modified Bessel function, 
A is the penetration depth, and fy = (r* — Tj)/\vi — r,-|. 
A cut off is placed on Ki(r/X) after it reaches an ex- 
tremely small value at r = 6A. The pinning force is 



f- P = E&C/pAp)^- r«"| 6((r p - | ri -r^|)/A)f 

(v) 

Here, O is the step function, rjj ' is the location of pin- 
ning site k, f p is the maximum pinning force, and = 
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The Lorentz force is modeled as a 



uniform force f^. All lengths, fields and forces are given 



in units of A, ^o/A 2 , and f 



Hi 



"A , respectively. 



We focus on experimentally accessible parameters that 
are close to those used in recent experiments 0. Fixing 
the sample size at 36A x 36A, with periodic boundary 



conditions, and the pinning radius at 



0.3A, we ex- 



amine an 18 x 18 square pinning array with N p = 324, 
giving a pinning density of n p = 0.25/A 2 . We slowly in- 
crease the driving force fd along the horizontal symmetry 
axis (x-axis) of the pinning lattice and compute the aver- 
age velocity in the x-direction, V x = (1/N V ) J2i=i v « ' 
This quantity is proportional to a macroscopically mea- 
sured voltage-current V(I) curve. In order to separate 
the different effects that each of the pinning and system 
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parameters have on the vortex dynamics, we fix all the 
parameters and vary only one at a time. 
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FIG. 1. Average vortex velocity V x versus driving force 
f d , for B/B^ = 1.062, f p = 0.625/ , r p = 0.3A, and 
B$ — 0.25$o/A 2 , with the pinning sites located in a square 
array; fd is increased from to 0.8/o. Several remarkable 
jumps in the curve can be clearly seen which correspond to 
transitions in the dynamical behavior of the driven lattice. To 
better identify the phases we have numbered them I-IV. The 
inset shows the hysteresis curve as fd is increased to 0.7/o 
and then decreased (phases now shown with unbold Roman 
numerals) to zero. Some phase boundaries, II — III and III-IV, 
show hysteresis while others do not. 

In Fig. 1 we present a typical V(I) curve for B > B^ as 
the driving force is linearly increased from to 0.8/0 with 
f p = 0.625/o, B^ = 0.25$ /A 2 , and B/B^ = 1.062. The 
V(I) curve exhibits several remarkable features which 
clearly appear as discontinuous jumps and drops in V x . 
We label each of the features as regions I through V, and 
approximate the fraction of flux lines which are mobile 
at a specific driving force with: a = V x / fa- As fd is in- 
creased in region I, V x is zero, indicating that the vortex 
lattice is pinned. At fd = 0.146/o, the onset of region 

II is marked by a finite V x , caused by the depinning of 
interstitial vortices. Here, a = 0.06, confirming that only 
the interstitial vortices are mobile since the percentage of 
vortices above B^ is also (B — B^/B^ — 0.06. Region 

III begins at fd = 0.406/o, where a very sharp jump up 
in V x is seen, along with an increase in the number of 
mobile vortices to a = 0.44. The velocity fluctuations 
5V X are also much larger. At fd = 0.462/ , region IV 
appears with a sudden drop in V x , with a = 0.23 and a 
reduction in SV X . Finally, at fd = 0.612/o, just under 
the pinning force of each pinning site (f p = 0.625/o), the 
entire lattice becomes mobile with a = 1, and the system 
enters region V. 

To further characterize the dynamic phases, we have 
performed a number of hysteresis runs where the driving 
force is slowly increased and then decreased. In the inset 
of Fig. 1 we show a typical curve for the system with the 
same parameters as in Fig. 1. There is little hysteresis 
for the transition IV-V. However, there is a very strong 



hysteresis at the transitions II— III and III-IV that per- 
sists in larger system sizes and also for low values of the 
spatial disorder. The hysteresis and the sharp jumps sug- 
gest that the II— III and III-IV phase boundaries might 
be first order. Details will be presented elsewhere Q |. 

We use a series of snapshots of the vortex position and 
vortex flow paths for regions II through V to show ex- 
plicitly that the features in the V(I) curve correspond to 
different plastic flow phases. Figure 2(a) shows the vor- 
tex trajectories in region II of the current-voltage curve 
from Fig. 1. It is clear that in region II only the intersti- 
tial vortices are mobile while the commensurate vortices 
remain pinned. The motion is confined to ID channels 
between the rows of pinning sites due to the square sym- 
metry imposed by the pinned vortices. Such flow behav- 
ior for vortices in square pinning arrays has recently been 
experimentally observed l|. 

Figure 2(b) illustrates that the vortex trajectories dif- 
fer greatly from region III to II. The vortex lattice is now 
disordered, and the flow is no longer ID but consists of 
channels that wind in both the x and y directions. Pin-to- 
pin motion also appears. Unlike the motion in region II, 
where only the interstitial vortices move and vortices at 
the pinning sites remain pinned, all the vortices in region 
III take part in the motion, with any one vortex moving 
for a time and then being temporarily trapped [ fl2"| |. 

Another significant change in the vortex motion ap- 
pears when the system enters region IV, as seen in 
Fig. 2(c). The vortex trajectories become more ordered 
and return to an exclusively ID flow, with the mobile vor- 
tices moving along the pinning rows rather than between 
the rows as in region II. Only certain rows are mobile, 
and in these rows the additional vortices above B^ leave 
their positions between the pinning rows to create mov- 
ing incommensurate ID structures along the pinning row. 
An entire row does not move continuously, but instead 
a pulse appears in which only four vortices, near the in- 
commensurate segment of the vortex row, are mobile. As 
this pulse or discommensuration moves across the sam- 
ple, each vortex in the row is displaced by a single pinning 
lattice constant a. This disturbance is thus crossing the 
sample much more rapidly than the vortices themselves. 

The vortex trajectories for region V, in which the entire 
lattice is moving, are shown in Fig. 2(d). Some portions 
of the vortex lattice have a distorted triangular order, al- 
though the incommensurabilities from region IV are still 
present. The flow remains strictly ID and along the pin- 
ning rows as in region IV, except that now all the rows 
are mobile. Rows with an incommensurate number of 
vortices move faster than the commensurate rows. As fd 
is increased further, the density of incommensurabilities 
and the vortex lattice structure do not change, so the 
system is always undergoing plastic flow and a moving 
crystal is not formed. 

The onset of these different phases is described with 
force balance arguments that take into account the cou- 
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pling of the two different species of vortices, interstitial 
and (pinned) commensurate. In region II, while the com- 
mensurate vortices remain pinned, the interstitial vor- 
tices begin flowing at a well defined driving force and 
exert a force f c -i C on the pinned commensurate vortices. 
The total force on a commensurate vortex thus consists 
of the driving force f<j, pinning force f p , and forces f c _ c 
from commensurate and f c -i C interstitial vortices. The 
commensurate vortex will remain pinned as long the fol- 
lowing inequality holds: 



If J > \fd 



fc-c 



(1) 



Since here we are using a square pinning lattice, from 
symmetry we have f c _ c = 0. If there are no incommensu- 
rate vortices, f c _j c = 0, and commensurate vortices depin 
at |f(j| = |£p|. When interstitial vortices are present, the 
term f c _i C causes commensurate vortices to depin before 
|f(i| = |f p |. These vortices depin more vortices so that 
the number of mobile vortices increases and region III 
appears. As long as |f<j| < |f p |, not all the vortices will be 
mobile so that a < 1. For the parameters used in Fig. 1, 
the density of interstitial vortices is sufficiently low that 
they do not interact significantly with each other. In this 
case, we can solve Eq. (1) for the transition from region 
II to region III to give 
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With the parameters for the onset of region III, Eq. (2) 
gives f p = 0.621/o, which is in very good agreement with 
the value of f p = 0.625/o used in the simulation. 

The appearance of ID motion exactly along the rows 
of pins in region IV might seem counterintuitive since 
for B/B,), > 1 and fd = 0, when the vortices are not 
moving, an incommensurate vortex located along a pin- 
ning row is unstable to perturbations in the y-direction 
and will fall into the interstitial area between rows. For 
moving vortices, the situation is quite different since the 
vortices spend part of their time in the pinning sites. The 
pinning sites create a stabilizing force against perturba- 
tions in the transverse direction, confining the motion 
along the pinning rows. When the density of interstitial 
vortices is low, the onset of region IV occurs when the 
driving force is strong enough that interstitial vortices 
can depin commensurate vortices from a distance a/2 in 
the longitudinal direction. This distance, and especially 
the repulsion from the remaining pinned vortices, allow 
the interstitial vortex to move towards the just-vacated 
pin site. For sufficiently strong fd, it will remain moving 
along the pinned rows JnJ. The transition from region 
III to region IV should occur at 
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where r% = (a/v2 + r p )/A. With the parameters used 
in Fig. 1 at the onset of the ID incommensurate flow, 
Eq. (3) gives f p — 0.624/ , which is in very good agree- 
ment with the numerical value for f p shown in Fig. 1. 

To better characterize the flow behavior, we systemat- 
ically vary f p with the rest of the parameters fixed. The 
resulting phase diagram in Fig. 3(a) outlines the onset 
of the different dynamical phases. As f p is increased, 
region I saturates at a value of fd ~ O.I46/o. This oc- 
curs because, although the pinning force f p of the pin- 
ning sites is being increased, the vortex-vortex interac- 
tions which determine the interstitial pinning force are 
not changed. Region II only occurs when f p > 0.37/o, 
since for f p < 0.37/o, Eq. (2) cannot be satisfied even 
for very low fd, and as soon as the interstitial vortices 
move they start to depin commensurate vortices. The 
same argument applies for the onset of region IV, which 
extends to even lower f p values. The II— III phase bound- 
ary follows Eq. (2), which for high fd goes as f p ~ fd, in 
agreement with the phase diagram. Similarly, the III IV 
boundary follows Eq. (3), which is also linear for large 
fd- The onset of region V also goes as f p cx fd- 

Next we vary the commensurability, from B/B^ = 0.75 
to 1.7, producing the phase diagram in Fig. 3(b). Just 
above B/B^ = 1, the five phases of Fig. 1 are present. As 
B/Bff, is increased, the disordered-flow region III grows 
while the ordered- flow regions II and IV shrink. This 
is expected since an increase in B/B^ effectively intro- 
duces more disorder via the addition of more interstitial 
vortices. For B/B^ > 1.3, the flow becomes more dis- 
ordered, the II— III phase boundary becomes ill-defined, 
and the flow in region V is no longer composed of ID 
incommensurate flow along the pinning rows but has a 
number of vortices flowing between the pinning rows. At 
B/Bff, = 1, the commensurate case, we find only two 
phases: pinned and flowing, with the onset of flow occur- 
ring at fd ~ fp- The vortex flow for the commensurate 
case is elastic since, unlike the case for B jB^, > 1, there 
are no discommensurations that cause certain rows to 
move faster. For B / B^ < 1, where a number of vacancies 
appear in the vortex lattice, we observe a new vacancy 
flow phase, marked region VI. The depining force for the 
onset of vacancy motion is considerably higher than that 
for the onset of interstitial motion, in agreement with 
experiments [0. Flow in V is like in V, but now the 
faster-moving rows have vacancies. 

Figure 3(c) shows a phase diagram in which the po- 
sitions of the pinning sites are gradually disordered by 
randomly displacing them up to an amount Sr from the 
ordered pinning lattice position. In terms of the pinning 
lattice constant a, we consider the case Sr = a/2 to be 
a good approximation to a random pinning array. The 
disordered region III grows and dominates the phase di- 
agram for disorder greater then Sr = a/6, so that only 
three phases occur, in agreement with other simulations 
of random pinning arrays [0j4|J 
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Figure 3(d) shows the phase diagram for increasing 
both vortex n v and pin n p density, determined by B,p. 
Here, the ordered-flow phase IV grows while the disor- 
dered flow region III shrinks. Increasing B^, increases n v 
and thus effectively decreases f p . As seen in Fig. 3(a), 
this suppresses II and III, while favoring IV. Phase dia- 
grams for varying pinning radius r p , angle of drive, tem- 
perature, and triangular pinning arrays will be presented 
and discussed in detail elsewhere Q. 

In conclusion, we have demonstrated several novel dy- 
namical phases which are very distinct from those found 
in random pinning arrays. These phases are marked by 
pronounced changes in the V(I) curves and noise fluc- 
tuations which should be very accessible experimentally. 
We have directly related these features to the pronounced 
changes in the vortex lattice structure and flow pat- 
tern. These phases result from the coupling between two 
species of vortices. We have shown that certain phases 
show strong hysteresis while others do not. The depen- 
dence of these phases on various pinning and system pa- 
rameters has been extensively studied and summarized 
in a series of phase diagrams [[yj. We hope that these 
results encourage the experimental search of these new 
dynamic phases. 
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FIG. 2. Trajectories of the east-bound flowing vortices, for 
the voltage-current curve in Fig. 1, for regions: II, interstitial 
flow (a); III, disordered flow (b); IV, incommensurate ID flow 
(c); and V, moving (in)commensurate rows (d). The vortices 
are represented by black dots and pinning sites by open cir- 
cles. For clarity, only a 20 A x 20 A subset of the 36 A x 36 A 
sample is shown. 
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FIG. 3. Dynamic phase diagrams with a square pinning 
lattice. Here, B/B^ = 1.062, B$ = 0.25$ /A 2 , r v = 0.3A, and 
f p — 0.625/o, unless otherwise noted, (a) Pinning force f v 
versus driving force fd- As fd is increased, the phase bound- 
aries II— III, III-IV, and IV-V become linear, (b) B jB^ versus 
fd- For B/Bj, > 1, regions I through V can be observed, with 
the disordered region III growing and the ordered-flow regions 
II and IV reducing in size. A similar situation occurs in (c) 
for gradually increasing the amount of disorder in the location 
of the pins, (d) Dynamic phase diagram for B^ versus fd- 
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